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Abstract 

In many practical situations, incertitudes affect the mechanical behaviour that is given 
by a family of graphs instead of a single one. In this paper, we show how the bipotential 
method is able to capture such blurred constitutive laws, using bipotential convex covers. 

1 Introduction 

The constitutive laws of the materials can be represented, as in Elasticity, by a univalued 
mapping T : X ^ y or, as in in Plasticity, can be put in the form of a multivalued operator 
T : X — > 2^. Equivalently, a constitutive law can be seen as the graph M C X x Y of the 
operator T, that is M = Graph{T) = {{x, y) ■ y & T{x)}. X and Y are spaces of dual variables, 
for example X may be a space of stresses and Y may be a space of deformation rates. The 
duality between these spaces is a function (■, ) : X x y ^ M. 

If the graph M is maximal cyclically monotone, then there is a convex and lower semi- 
continuous (l.s.c.) function (j) : X ^ M., called a superpotential (or pseudo-potential), such 
that M — Graph{d(f)) , where d(j) is the subdifferential of (j) and the function cj) is determined by 
the graph M, up to an additive constant. 

The constitutive laws of dissipative materials admitting a superpotential can be put into the 
form y G d<j){x). Such constitutive laws are often qualified as standard [9] and the law is said 
to be a normality law, a subnormality law or an associated law. However, many experimental 
laws proposed these last decades, particularly in Plasticity, are non associated. For such laws, 
we proposed in [13j a suitable modelization thanks to a function called bipotential. 

The laws admitting a bipotential are called laws of implicit standard materials because they 
have the form y G y){x) , which is a subnormality law but the relation between x and y is 
implicit. 

The bipotential theory allows, in connection with the calculus of variation, to model a wide 
spectrum of non associated constitutive laws. Examples of such non associated constitutive laws 
are: non-associated Driicker-Prager [15] and Cam-Clay models 16 in soil mechanics, cyclic 
Plasticity ([14],[2]) and Viscoplasticity [W of metals with non linear kincmatical hardening 
rule, Lemaitre's damage law [Ij, the coaxial laws ([l7j,[19j), the Coulomb's friction law [13], 
[H], [3], [I], [8], [11], [H], [H], [ig. a complete survey can be fomid in [17 . In the previous 
works, robust numerical algorithms were proposed to solve structural mechanics problems. 
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The cornerstone inequality in the definition of bipotentials (definition 12.11 (b)) extends the 
Fenchel's inequahty. In particular, to any superpotential (p is associated the separable bipo- 
tential: 

6(x,y) = 0(x)+r(y) (1.0.1) 

where (/)* is the Fenchel conjugate of (}> (with respect to the duality between the spaces X 
and Y). The implicit subnormality law y G db{-,y){x) becomes the associated law y € d(j){x). 
However, there are many bipotentials which can not be expressed in the form (jl.O.ip . 

For all the particular constitutive laws previously mentioned, the bipotentials were heuristi- 
cally constructed, without knowing beforehand the conditions under which ones the law admits 
a bipotential, nor a systematic algorithm to construct this bipotential. 

In [4] we solved two key problems: (a) when the graph of a given multivalued operator can 
be expressed as the set of critical points of a bipotentials, and (b) a method of construction of 
a bipotential associated (in the sense of point (a)) to a multivalued, typically non monotone, 
operator. The main tool was the notion of convex lagrangian cover of the graph of the 
multivalued operator, and a related notion of implicit convexity of this cover. The results of [3] 
apply only to bi-convex, bi-closed graphs (for short BB-graphs) admitting at least one convex 
lagrangian cover by maximal cyclically monotone graphs. This is a rather large class 
of graph of multivalued operators but important applications to the mechanics, such as the 
bipotential associated to contact with friction |13| . are not in this class. 

In more recent papers [5] , [5] , we proposed an extension of the method presented in [1] to 
a more general class of BB-graphs. This is done in two steps. In the first step we proved that 
the intersection of two maximal cyclically monotone graphs is the critical set of a bipotential 
if and only if a condition formulated in terms of the inf convolution of a family of convex Isc 
functions is true [6] . In the second step we extended the main result of [14] by replacing the 
notion of convex lagrangian cover with the one of bipotential convex cover (definition 13. 2p . 
In this way we were able to apply our results to the bipotential for the Coulomb's friction law. 

The purpose of this paper is to describe a new application of bipotentials. In many practical 
situations, incertitudes affect the mechanical behaviour. In other words, we tolerate indeter- 
minacy of the constitutive law which is represented by a family of graphs instead of a single 
one. In particular, when no solution can be found for ill-posed problems, relaxation of stronger 
conditions on the material behavior would allow to provide at least an approximate solution. 
Our aim now is to show how the bipotential is able to capture such blurred constitutive law, 
by using bipotential convex covers. 

The main results of the paper are Propositions l4.HI5.1l andl 6.1l In Proposition |4T] we find a 
bipotential for a blurred Elasticity law, or equivalently, we show that such law can be expressed 
as a implicit subnormality law. We pass then to a more difficult blurred Plasticity law, with 
a variable yielding threshold rj taking arbitrary values in [A_, A+]. For this law we achieve the 
same as previously in Proposition l5.ll 

The third result concerns a blurred Coulomb friction law. The law of unilateral contact with 
Coulomb's dry friction is a typical example of a non associated constitutive law in mechanics 
which admits a bipotential [13j . It is to be remarked that the Coulomb friction law does not 
have a separated bipotential, because the graph of the law is not even monotone. In Proposition 
16.11 we are able to provide a bipotential formulation for a blurred Coulomb friction law with 
arbitrary values for the friction coefficient /i in a range ^+]. 

Aknowledgements. The first two authors acknowledge the support from the European Asso- 
ciated Laboratory " Math Mode" associating the Laboratoire de Mathematiques de I'Universite 
Paris-Sud (UMR 8628) and the "Simion Stoilow" Institute of Mathematics of the Romanian 
Academy. The first author thanks to M. Jean and M. Raous of the Laboratoire de Mecanique 
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and the second one for the kind invitation to give a seminar which led to discussions where this 
idea was suggested. 



2 Bipotentials 

X and Y are topological, locally convex, real vector spaces of dual variables x G X and y ^Y, 
with the duality product (•,•) : X x Y ^ R. We shall suppose that X,Y have topologies 
compatible with the duality product, that is: any continuous linear functional on X (resp. Y) 
has the form x ^ {x,y), for some y € Y (resp. y i— s- {x,y), for some x G X). We use the 
notation: R = M U {+00}. For any convex and closed set A C X, its indicator function, xa, is 
defined by 

. . f if X e A 
Xa[x) <^ otherwise 

The subgradient of a function : X ^ M at a point x ^ X is the (possibly empty) set: 

d4>{x) ^{ueY \ yzeX {z~x,u) < (j){z) - (j){x)} . 

Definition 2.1 A bipotential is a function b : X x Y ^ M., with the properties: 

(a) b is convex and lower semicontinuous in each argument; 

(b) for any x ^ X,y ^Y we have b(x,y) > {x,y); 

(c) for any [x^y) ^ X xY we have the equivalences: 

yedb{;y){x) ^ xedbix,-){y) ^ b{x,y) = {x,y) . (2.0.1) 

The graph of b is 

M{b) = {{x,y) eXxY \ b{x,y) = {x,y)} . (2.0.2) 



If the graph M of a law is the graph of a bipotential 6, we say that the law (the graph) 
admits a bipotential. 

For any graph M <Z X xY , we can introduce the sections M{x) — {y &Y \ {x,y) € M} 
and M*{y) = {a; S X | (x, y) € M}. Hence the operator T assigns to each x E X the section 
M{x) and the inverse law assigns to each y G Y the section M*{y). 

Let a constitutive law be given by a graph M . Does it admit a bipotential? The existence 
problem is easily settled by the following result. 

Theorem 2.2 Given a non empty set M G X xY , there is a bipotential b such that M = M(b) 
if and only if for any x € X and y GY the sections M{x) and M*(y) are convex and closed. 

The proof can be found in 4 . Then we say that M is bi-convex and bi-closed, or in short 
that M is a BB-graph. 

If the law is represented by a BB-graph, then a closely related topic is to know whether 
the bipotential is unique. The answer is no. The proof of the previous result is based on the 
introduction of the bipotential 

boo{x,y) = {x,y) +XM{x,y). 

If the graph M = Graph(d(j)) is cyclically monotone maximal, then it admits at least two 
distinct bipotentials: the separable bipotential defined by 11.0. II and boa ■ Therefore the graph 
of the law alone is not sufficient to uniquely define the bipotential. 
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3 Bipotential convex covers 



Theorem 12.21 does not give a satisfying bipotential for a given multivalued constitutive law, 
because the bipotential boo is somehow degenerate. We would like to find a method of con- 
struction of bipotentials which for a given BB-graph M will return a bipotential b which is not 
everywhere infinite outside the graph M , and such that if M is maximal cyclically monotone 
then the method will give us a separable bipotential. 

We saw that the graph alone is not sufficient to construct interesting bipotentials. We need 
more information to start from. This is provided by the notion of bipotential convex cover. 

Let Bp{X, Y) be the set of all bipotentials b : X x Y ^ R. We shall need the following 
definitions (4.1 and 4.1 [6]). 

Definition 3.1 Let A be an arbitrary non empty set and V a real vector space. The function 
f : A X V ^ M. is implicitly convex if for any two elements (Ai, zi), (A2, Z2) E A x V and for 

any two numbers a, (3 E [0, 1] with a + P = 1 there exists A G A such that 

f{\azi+[iz2) < af{XuZi) + l3f{X2,Z2) . (3.0.1) 



Definition 3.2 A bipotential convex cover of the non empty set AI is a function A G A 
b\ from A with values in the set Bp{X,Y), with the properties: 

(a) The set A is a non empty compact topological space, 

(b) Let / :AxXxy^MU {+cx3} be the function defined by 

f{\x,y) = bx{x,y). 

Then for any x E X and for any y E Y the functions f{-,x,-) : A x Y ^ M. and 
f{-,-,y) : A X X — > R are lower semi continuous on the product spaces A xY and 
respectively A x X endowed with the standard topology, 

(c) We have AI = [j M{bx). 

AgA 

(d) with the notations from point (h), the functions /(•, x, ■) and /(•, •, y) are implicitly convex 
in the sense of Definition \3.1\ . 

The next theorem, [6] theorem 4.6, is the key result needed further. 

Theorem 3.3 Let X 1-^ b\ be a bipotential convex cover of the graph M and b : X x Y —> R 
defined by 

b{x,y) = ini{bx{x,y) \ X E A} . (3.0.2) 
Then b is a bipotential and M = M{b). 

The result is rather surprising because an inferior envelop of functions, even convex, is not 
generally a convex function. The property (d) of the Definition 13.21 is essential to ensure the 
convexity properties of b. 
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4 Application to a Elasticity law with thick line 



The finite element method is a numerical method of discretizing continua, widely used today 
for solving structural mechanics problems. The accuracy of the approximative solution can be 
controlled by computing a posteriori estimators. Three different approaches were proposed. 
Contrary to the methods based on the equilibrium residuals ([10], [H], [52], [53], [15 ), the 
ones using smoothing techniques ([IS], [5S], [57]) and the dual analysis based on upper and 
lower bounds for the energy ([28], [29]), the method of the error on the constitutive law (or 
constitutive relation error) is based on mechanical concepts and can be more naturally extended 
to non linear problems of evolution ([30], [31], [32], [33], |34], [35], |36], [37]). We use such an 
idea in the sense that we admit an error on the elastic law instead of satisfying it exactly. 
We shall use bipotentials constructed from bipotential convex covers, in order to formulate the 
elasticity law with thick line as a implicit subnormality law. 

We shall take X = Y = M" and the duality product is the usual scalar product in M". 
let us consider the elastic linear law y = Xx with A > which is the most simple example 
of linear elastic law where the dual variables x and y are vectors. In the present application, 
the material parameter, the "elastic modulus" A, has a fixed value but we allow some error on 
the constitutive law a = y — Xx with a fixed margin of tolerance e > 0. In other words, we 
considered the blurred elastic law described by the BB-graph: 

M = {{x,y) £ X xY \ \\y ~ Xx \\<e} . (4.0.1) 




(a) (b) 
Figure 4.1: Elasticity law. (a): ideal law with thin line, (b): law with thick line 

In figure 14. 1[ a pictural representation of the blurred law is a graph with " thick line" (dis- 
played at the right) in contrast to the ideal law of the previous section with a "thin line" 
(displayed at the left). 

Proposition 4.1 The blurred elastic law with the graph M J[ ) is represented by the bipo- 
tential: 

bix,y)^{x,y) + j^ {(\\y-Xx\\~e)^)' (4.0.2) 

with the notation (a)^ = max(Q!, 0). 

Equivalently, the relation \\y — Xx\\ < e can be put in the form y G db(-,y)(x), with b the 
bipotential defined by ^J^UT^. 
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Proof. The graph M is clearly bi-convex and bi-closed. We construct a bipotential convex 
cover of M by assigning to the parameter a e i?(e) the set 

Ma = {{x,y) e X X r I y^\x} , 

which can be seen as the graph of the elastic law y = Xx + a with an "initial stress" a. it is 
clear that Ma is the subdifferential d(f>a of the potential: 

(jJaix) = II a; f + {x,a) . 



Its Fenchel conjugate is: 



0:(y)-^ll2/-af 



Let {ba : a £ B{e)} be the collection of the separated bipotentials: 

ba{x,y) =Mx) + ra{y) = ^W^W^ +{x,a) + jj-\\y-a\W 

ba{x,y) = {x,y) + ^\\y-a-Xx\\'^ 

We want to verify that it defines a bipotential convex cover of M. The conditions (a) to (c) of 
Definition 13.21 are fulfilled. We have to prove the last condition (d). For the implicit convexity 
of f{-,x,-), the incqualitv 13.0.11 is then: for any x £ X, yi,y2 G Y, ai,a2 G B{e), a, (3 G [0,1], 
a + /? = 1, there exists a E B{£) such that: 

a\\yi-ai - Xx |p +/3 || t/2 - a2 - Ax |p>|| ayi + (3y2 - a - Xx |p . 

Because the square of the norm is convex, an obvious choice for a is a = aai + (3a2 ■ For the 
implicit convexity of /(■,■,?/), the demonstration is similar. We apply Theorem 13.31 and we 
obtain a bipotential of the graph M with the form: 

b{x,y) = ini{ha{x,y) : a £ B{e)} = mi \^ \\ x f + {x,a) + ^ \\ y ~ a f : ||a||<e) 



^2 ' ' ' 2A 
We shall prove now that h has the form (|4.0.2p . Two events have to be considered: 

- The infimum is realized at a such that || a \\< e. Hence ba{x, y) is stationary with respect 
to a, that is: a — y — Xx. Eliminating a by this relation, a straightforward calculation 
shows that the infimum is: 

b{x,y) = {x,y) . 

- Otherwise, introducing a Lagrange multiplier ry, we have: 

b{x,y) = sup{inf{ baix,y) + r/ {\\ a f -e^) : a G F} : t] > O} 

The stationarity condition with respect to a gives: 

y- Xx 



1 + 27/ ■ 

Hence the constraint || a \\ — e allows to deduce the value of the Lagrange multiplier: 



(4.0.3) 



v = U-\\y^^x\\-i] ■ (4.0.4) 
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Introducing expression 14 . . 41 into 14 . . 3l leads to: a = e-r-^ — rr. Eliminating a by this relation 

II 2/ -Ax II 

gives the value of the infimum: 

1 2 

b{x,y) ^ {x,y) + — {\\y- Xx\\-e) . 

We verify immediately that b has indeed the expression ()4.0.2I □ 

The same reasoning may be performed for the more realistic elasticity law y — Kx with K 
a elasticity tensor, but the computations are more involved. 



5 Application to a blurred Plasticity law 

We want now to extend the previous ideas to non smooth constitutive laws. Let us consider 
the Plasticity law with a yielding threshold A for which the plastic domain is the the closed ball 
B{X) of center and radius A: B{\) = {y E Y : \\y\\ < A}. The Plasticity law is given by the 
maximal cyclically monotone graph: 

^ {iO,y)€X xY : \\y\\<X} U {{x,y)eXxY: \\ y \\= \,3r, > 0,x ^ r,y} . 

The graph M\ is the subdifferential d4>\ of the potential: 

0a(x) = A||x|| . (5.0.1) 

Its Fenchel conjugate is: (f>X{y) — Xb{\){v)- We allow some error on the constitutive law: 

y e dcj)x(x) + a 

with a fixed margin of tolerance e > on the norm of a. In other words, we consider the blurred 
plastic law described by the graph: 

M = {{x, y)eXxY I 3ae B{t), y £ d(l)x{x) + a } . (5.0.2) 



X 



■x_ 



(a) (b) 
Figure 5.1: Plasticity law. (a): ideal law with thin line, (b): law with thick line 

In figure 15. 1[ a pictural representation of the blurred law is a graph with " thick line" (dis- 
played at the right) in contrast to the ideal law with a "thin line" (displayed at the left). 
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Proposition 5.1 The blurred plastic law with the graph M i f 5. 0. 2\) is represented by the bipo- 
tential: 

b{x,y) = sup{X-, II y II) II a; || +XBix^){v) (5.0.3) 

with the notation A± = A ± e. 

Equivalently, the relation y G d(j){x) + a, \\a\\ < e, can be put in the form y G db{-,y){x), 
with b the bipotential defined by i5. 0.3]) . 

Proof. As previously, we intend to construct a bipotential for M from a bipotential convex 
cover. This time a is not the good parameter for a bipotential convex cover because it is difficult 
to check all the conditions of Theorem 13.31 for this parameter. 

The bipotential convex cover will have as parameter a variable yielding threshold 77 taking 
arbitrary values in [A_, A+] where A± = A ± e. The graph M (|5.0.2I ) admits the description: 

A/ = {(0,2/) e X X r I II y ||< A_}U (5.0.4) 

U{{x,y) eX xY I A- <|| y ||< A+ and 3r]>0, x = rjy} 

Let {bjj : 77 e [A_, A+]} be the collection of separated bipotentials: 

f{v,x,y)^brj{x,y) = ??||a;|| + Xi3(,,) (y) ■ 

We shall verify that this is a bipotential convex cover of M. Indeed, the conditions (a) to (c) 
of Definition 13.21 are obviously fulfilled. With f{r],x,y) = bjj{x,y), we have to prove the last 
condition (d). For the implicit convexity of /(•, •, y) we have to prove that for any xi,X2 G X , 
y ^Y, T]i,ri2 e [A_, A+], a, /? e [0, 1], a + /3 = 1 there exists rj G [A_, A+] such that: 

aT]i II 0:1 II +Pl]2 II X2 II +XB{m)iy) +XS(r,2)(2/) > V II "2:1 +1^X2 II +XB{r^){y) ■ 

We choose rj — min(?7i, 772). The previous inequality becomes: for any || y ||< 77 

ar^i II x\ II +/37/2 II X2 ||> r\ || oix\ + (3x2 \\ , 

which is true because of the convexity of the norm. 

For the implicit convexity of f{-,x,-) we have to prove that: for any x G X, yi,y2 G Y, 
ryi, 772 G [A_, A+], a, /3 G [0, 1], a + /? = Ithere exists rj G [A_ , A+] such that: 

{a'ni+f3r]2) II X II +XB{rn)iyi) + XB{r,2){y2) > ?7 II a; II +XB(r,){ayi + f3y2) 
This time we choose — arji + j3ri2. Indeed, for any || yi ||< Ai and || y2 ||< A2, we have: 

II ayi + /3y2 II < a II yi || +/3 || y2 || < am + PV2^V 
therefore the inequality we want to prove becomes trivial: 

{aril + (3ri2) \\ x \\> f] \\ x \\ . 
Hence we are in the conditions of applying Theorem l3.3l The graph M admits the bipotential 
b{x,y) ^ inf {A||a;|| +Xi3(A)(y) : Ag[A_,A+]} . 
In order to compute 6, three events have to be considered: 

- for II y \\> K, b{x,y) = +00 , 

- for A_ <|| y ||< A+, b{x,y) = mi {X\\x\\ : || y ||< A < A+} - ||y||||a:|| , 

- for II y ||< A_, b{x,y) = inf {A||a;|| : A_ < A < A+} = A_||a:|| . 
The proof is done. □ 
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6 Application to a blurred Coulomb's friction law 



The law of unilateral contact with Coulomb's dry friction is a typical example of what is called 
a non associated constitutive law in mechanics. Despite of its rather complex structure, it is 
worthwhile to have interest in it because of its importance in many practical problems. 

We shall not discuss here the phenomenal and experimental aspects but only the mathe- 
matical modeling with respect to the bipotential theory. To be short, the space X = M'^ is the 
one of relative velocities between points of two bodies, and the space Y, identified also to R'^, 
is the one of the contact reaction stresses. The duality product is the usual scalar product. We 
put 

{xn,xt) e X ^RxR\ {yn,yt) eY = RxR\ 

where x„ is the gap velocity, xt is the sliding velocity, z/n is the contact pressure and yt is minus 
the friction stress. The friction coefficient is /i > 0. The graph of the law of unilateral contact 
with Coulomb's dry friction is defined as the union of three sets, respectively corresponding to 
the 'body separation', the 'sticking' and the 'sliding'. 

M^^{ix,0)eX xY I Xn<0}U{iO,y)eXxY \ \\ yt \\< f^Vn} U (6.0.1) 



\J<{x,y) £ X xY I Xn=0, xt^ 0, yt = m 




It is well known that this graph is not monotone, then not cyclically monotone. As usual, we 
introduce Coulomb's cone 



Kf, = {{yn,yt) eY \ \\yt\\< f^yn} , 

and its conjugate cone 

= {{Xn,Xt) £ X I Xt \\ +Xn < 0}. 

In particular, we have 

Ko = {(y„,0) e F I 2/„ > 0} , K* = {{x^,xt) eX \ x,, < 0} . 

In [5], [6], we obtained as application of Theorem 13 . 31 the following expression: 

b,,{x,y) = fiyn II Xt II +XK^{y) + XK'ix) (6.0.2) 

recovering the bipotential firstly obtained in heuristic way in |13| . 

We shall modify the unilateral contact law with Coulomb's dry friction by allowing arbitrary 
values for the friction coefficient fi in a, range that leads to the blurred friction law 

represented by the graph: 

M = {{x,0) e X xY I x„<0}U (6.0.3) 
U{{0,y) e X xY I ^J.^yn <\\ yt \\< ^J.+yn, 3?7 > O, xt=riyt} 

Proposition 6.1 The blurred plastic law with the yraph M !i6.0.3\} is represented by the bipo- 
tential: 

b{x,y) = sup(/i_t/„, II yt \\) \\ xt \\ +XK^^iy) + Xk*{x) (6.0.4) 

Equivalently, the relation y G dh^{-,y), fi G with b^ defined by \6. 0.S\) . can be put in 

the form y e db{-,y){x), with b the bipotential defined by |g. 0.4^ . 
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Proof. The collection {6^ : fi G of bipotentials 6^ defined as in (|6.0.2p is a bipoten- 

tial convex cover. The demonstration of the conditions of Definition 13.21 is similar to the one of 
the blurred plasticity law and is not reproduced here. 

Applying Theorem 13.31 the graph (j6.0.3|) of the blurred law admits the bipotential: 

b{x,y) = inf {/ij/„ || Xt \\ +XK^{y) ■ fJ- e [n-,fi+]} + Xk*{x) . 

We shall compute this bipotential. If a; G Kq , three events have to be considered: 

- for II yt II > fi+Vn, b{x,y) = +oo , 

- for ^_y„ <|| yt \\< fi+yn, b{x,y) = inf {pi?/„||xt|| : || yt \\< fiyn < tJ.+yn} = btlllktll , 

- ioi \\ yt \\< fi^yn, b{x,y) = inf {^y„||a;t|| : /i_ < ^ < = ^_y„||xt|| . 

If X does not belongs to Kq , b{x, y) = +cxd . In short, the blurred Coulomb friction contact 
law given by the graph 16.0.31 admits the bipotential (|6.0.4p . □ 

7 Conclusion 

We showed how the bipotential method is able to capture blurred constitutive laws. The 
blurring of a constitutive law enters by allowing some error in the law with a fixed margin 
of tolerance. For example a material parameter takes indeterminate values in a compact set 
modeling some experimental incertitudes. The notion of bipotential convex cover leads to the 
construction of a bipotential for such a blurred constitutive law. These bipotentials could be 
used to represent structural mechanics problems with blurred constitutive laws by means of 
variational inequalities. 
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